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ABSTRACT 

We have conducted a systematic survey for intervening O VI absorbers in available echelle spectra 
of 16 QSOs at ZQso = 0.17 - 0.57. These spectra were obtained using HST/STIS with the E140M 
grating. Our search uncovered a total of 27 foreground O VI absorbers with rest-frame absorption 
equivalent width VFr.(1031) > 25 mA. Ten of these QSOs exhibit strong OVI absorbers in their 
vicinity. Our OVI survey does not require the known presence of Lya, and the echelle resolution 
allows us to identify the O VI absorption doublet based on their common line centroid and known flux 
ratio. We estimate the total redshift survey path, A z, using a series of Monte-Carlo simulations, and 
find that Az — 1.66,2.18, and 2.42 for absorbers of strength Wr — 30,50 and 80 mA, respectively, 
leading to a number density of dJ\f{W > 50 uik)/dz = 6.7 ± 1.7 and dN'CW > 30 mk)/dz = 
10.4 ± 2.2. In contrast, we also measure dAf/dz = 27 ± 9 for OVI absorbers of Wr > 50 mA at 
|Au| < 5000 kms~ from the background QSOs. Using the random sample of OVI absorbers with 
well characterized survey completeness, we estimate a mean cosmological mass density of the O VI 
gas r2(0^+) h = (1.7 ± 0.3) x 10~^. In addition, we show that < 5% of OVI absorbers originate in 
underdense regions that do not show a significant trace of H I. Furthermore, we show that the neutral 
gas column A^(II I) associated with these O VI absorbers spans nearly five orders of magnitude, and 
shows moderate correlation with A^(0 VI). Finally, while the number density of OVI absorbers varies 
substantially from one sightline to another, it also appears to be inversely correlated with the number 
density of H I absorbers along individual lines of sight. 

Subject headings: cosmology: observations — intergalactic medium — quasars: absorption lines 



1. INTRODUCTION 

Uncovering the missing baryons in the present-day uni- 
verse is one of the central issues in observational cosmol- 
ogy, because an accurate accounting of the baryon bud- 
get at different locations sets an important constraint 
in modeling galaxy formation and evolution. Over the 
past few years, measurements of the baryonic mass den- 
sity from independent methods have approached a con- 
sistent value with high precision, Ui, = 0.020 ± 0.002 
(iBurles et alTgOOiySt^el et al.lf2003l : lWang et al.ll2003l : 



lO'Meara et al.ll2006f ). But while > 95% of all baryons 
can be acc ounted for in the Lya forest at redshift z ~ 3 
(| Ranch et al. 1997), only 1/3 have been identified at red- 
shift z = in known components, such as stars and neu- 
tral gas in galaxies or hot plasma in galaxy clusters (e.g. 
iPersic fc Saluccilll992l : fFukugita et al.lll998f) . 

Cosmological simulations indicate that « 30 — 40% of 
all baryons at the present epoch reside in a diffuse warm- 
hot intergalactic medium (WHIM), with temperatures in 
the range 10*^ < T < 10^ K, rather than in virialized ha- 
los where galaxies or g alaxy groups lie (jCen fc Ostriked 
Il999t iDave et aLll2001[ ). The simulations also show that 
this WHIM is shock-heated as tenuous gas accretes onto 
large-scale filamentary structures and remains hot owing 
to the low gas density and inefficient cooling. Locat- 
ing the WHIM in order to verify these models is diffi- 
cult, because the temperature of the gas is too low to 

^ Based in part on observations with the NASA/ESA Hubble 
Space Telescope, obtained at the Space Telescope Science Institute, 
which is operated by the Association of Universities for Research 
in Astronomy, Inc., under NASA contract NAS5-26555. 



be easily detected with the current generation of x-ray 
satellites and is too hot to be observed from the ground 
using optical-IR facilities. Simple ionization modeling 
of hot gas at T ^ 3 X 10^ K and metallicity 0.1 solar 
has shown that the gaseous clouds would imprint promi- 
nent absorption features produced by ionic transitions 
such as OVI AA 1031, 1037 and Ne VIH A A 770, 780 in 
the spectra of backgro und QSOs (e.g. lVerner et al.lll994l : 
iMulchaev et al.lll996l ). Surveys for the WHIM based on 
the presence of these ions at 2 < 1 therefore require a 
h igh-resolutio n UV s pectrograph in space. 

iTripp et"an ()2000l ) searched for intervening OVI ab- 
sorbers along the sightHne toward H1821-I-643 (zqso — 
0.297) in echelle spectra {R = A/AA « 43,000) ob- 
tained using the Space Telescope Imaging Spectrograph 
(STIS) on board the Hubble Space Telescope (HST). 
The search indicated that intervening O VI absorp- 
tion systems found in QSO spectra may indeed con- 
tain a significant fraction of baryons that had been 
missed by surveys based on the direct detection of pho- 
tons emitted by known objects. Subsequent searches 
along individual sightlines have b een carried o ut to- 
ward PG0953-f415 at zqso = 0.23 9 (ISavaee et al.l [2002'). 



PG1259+593 at zqso = 0.478 llRicht er et aLr |2004): 
PG1116-h215 at ZQSO = 0.1765 (Semba ch et "all 20041. 
PKS0405~123 at zqsq = 0.573 (Prochaska ct al.' 2004), 
and HE0226-4110 at zqso = 0.495 (Lchner et al. 2006). 
The measured line density of intervening O VI absorbers 
with rest-frame absorption equivalent width of the 1031 
transition Wr > 50 mA ranges from dM{OV\) / dz = 11 
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to 36^ per Une of sight at 2:0 vi = [0,0.5]. 

The large scatter in the observed number density of 
O VI absorbers suggests a significant variation in the spa- 
tial distribution of 0^+ ions across different sightlines. It 
also underscores a potential systematic bias in evaluat- 
ing the cosmological mass density of warm-hot gas based 
on selective quasar sightlines. The first blind search of 
OVI absorbers was carried out at zpyi = [0.5, 2 1] to - 
ward 11 background quasars by iBurles fc Tvtlerl (jl996l ) 
using medium resolution {R « 1300) and S/N « 20 
Faint Object Spectrograph spectra. They reported 
dAf{OVl)/dz = 1.0 ± 0.6 per line of sight at a mean 
redshift (zqvt) = 0.9 w i th W r > 210 ml. More re- 
cently, iDanforth &: Shull (j2005f ) conducted a search of 
O VI absorption features in known Lya absorbers at 
Zahs < 0.15, using Far Ultraviolet Spectroscopic Explorer 
(FUSE) spectra {R w 15, 000) of 31 active galactic nuclei 
(AGN). Their search yielded dA/'(0 VI)/dz = 9 ± 2 and 
17 ± 3 per hue of sight with Wr > 50 and Wr > 30 mA, 
respectively, at zoyi < 0.15, which is lower than mea- 
surements from individual lines of sight discussed earlier. 

Despite the large scatter in the number density of inter- 
vening O VI absorbers, these measurements imply that 
at least ^ 10% of the total baryons reside in the gaseous 
clouds probed by the O VI transition. Assuming a mean 
metallicity of 1/10 solar and an ionization fraction of 20% 
for the O VI absorbers, Tripp et al. (2000) derived a lower 
limit to the cosmological mass density in the ionized gas 
of r2b(0VI) > 0.003/1-1 Danforth & Shull (2005) 
derived r2f,(0VI) « 0.002 h^^. This general inference of 
O VI absorbers being a substantial reservoir of missing 
baryons is supported by numerical simulations that re- 
sult in a mean metallicity of 1/10 sol ar for the WHIM 
(ICen et al.l 120011 : iFang fc BrvanI 1200 iD . However, addi- 
tional uncertainties arise due to the unknown ionization 
state and metallicity of the gas. 

We have conducted a systematic survey of intervening 
OVI absorbers in available echelle spectra obtained us- 
ing HST/STIS with the E140M grating. A search of the 
HST data archive yielded high-resolution (i? ss 45, 000) 
spectra of 16 distant quasars at ^qso = 0.17 — 0.57 
with a mean S/N > 7 per resolution element. These 
quasars have been selected by different authors for dif- 
ferent studies and have bright UV fluxes that are suit- 
able for STIS echelle observations. Therefore, they rep- 
resent a random, flux-limited sample of quasar spectra 
for an unbiased search of intervening absorbers. The 
entire STIS echelle spectroscopic sample of quasars sum- 
marized in Table 1 defines the random sightlines along 
which our blind search of intervening O VI absorbers was 
conducted. 

The primary objectives of the blind OVI survey are: 
(1) to obtain a statistically representative sample of O VI 
absorbers along random lines of sight for an accurate 
measurement of the cosmological mass density of warm- 
hot gas at low redshift; and (2) to study the incidence 
of OVI absorbers along individual sightlines for quanti- 
fying the spatial variation of O VI absorbing gas. These 

^ ITripp et all pOOOh reported (W{OVl)/dz ~ 48 based on four 
OVI absorbers found with Wr > 30 mA toward H1821+643. 
Three of the four absorbers have Wr > 50 mA. We therefore scale 
dAf{OVl)/dz accordingly and work out dAf{OVT)/dz = 36 for 
OVI absorbers stronger than Wr = 50 mA. 



OVI absorbers identified along random sightlines also 
form a uniform sample for studying the physical nature 
of tenuous gas probed by the presence of 0^+ and for in- 
vestigating their large-scale galactic environment. Here 
wc focus our discussion on the statistical properties of 
OVI absorbers. Studies of the gas properties and dis- 
cussions of individual absorbers are given in an accom- 
panying paper (Thom & Chen, 2008; hereafter paper II), 
and studies of the surrounding galaxies will be presented 
in future papers. 

This paper is organized as follows. In Section 2, we 
summarize the available quasar spectra and the data 
quality. In Section 3, we describe our approach for iden- 
tifying an O VI AA 1031, 1037 absorption doublet with no 
prior knowledge of the presence of H I or other ionic tran- 
sitions. A Monte-Carlo analysis was performed for eval- 
uating the completeness of the survey. In Section 4, we 
present a statistical sample of O VI absorbers identified 
along the sightlines toward the 16 background quasars. 
We measure different statistical quantities of the ionized 
gas. In Section 5, we compare our search results with oth- 
ers and discuss the implications in Section 6. We adopt 
a A cosmology, VIm = 0.3 and = 0.7, with a dimen- 
sionless Hubble constant h = ifo/(100 km s^^ Mpc^"'^) 
throughout the paper. 

2. THE STIS ECHELLE SPECTRA 

The goal of our project is to conduct a blind sur- 
vey of OVI AA 1031, 1037 absorption doublet features 
along random lines of sight toward background quasars 
for establishing a statistically representative sample of 
the OVI absorbers at z < 0.5. The survey therefore 
requires high-quality UV spectra of distant quasars. In 
addition, we focused our search on available echelle qual- 
ity spectra in order to uncover the majority of O VI 
absorbing gas. Both numerical simulations presented 
by Chen et al. (2003) and the recent FUSE survey by 
iDanforth fc Shulll (120051 ) show that the bulk (> 60%) 
of O VI absorbers have rest-frame absorption equivalent 
width of the 1031 transition = 30 - 100 niA. An 
echelle spectrum with resolution element 6 v w 7 kms^^ 
at moderate S/N (w 5) per resolution clement is there- 
fore necessary for detecting these weak lines. This crite- 
rion excludes Faint Object Spectrograph data (typically 
(5w « 230 — 270 kms^ ) from our search. Finally, the 
throughput of STIS declines steeply at A < 1180 A, cor- 
responding to a minimum zqvi ~ 0.14, and we there- 
fore include in our sample only those sightlines toward 
quasars at zqso 0.16 for the search of intervening O VI 
absorbers. 

Our search in the HST data archive yielded 16 quasars 
observed with the STIS E140M grating. A summary 
of the STIS observations is presented in Table 1, where 
we list in columns (l)-(6), the QSO field, zqsoi mini- 
mum and maximum redshifts over which the sample of 
random OVI absorbers is established (^min, ^max), total 
exposure time, and the program ID. We note that Zmax is 
defined for O VI absorbers at velocity separation > 5000 
kms-i from the background QSO, but the line search is 
conducted through the emission redshift of the QSO. In- 
dividual flux-calibrated spectra were pr ocessed using the 
latest calibration files (|Aloisi et al.|[2005il and retrieved 
from the MAST server. For each QSO, individual echelle 
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TABLE 1 

Summary of the STIS Echelle Spectra 



QSO 2QSO ^min Zmax'' toxp PID 

(1) (2) (3) (4) (5) (6) 



HE0226-4110 0.4950 0.1154 0.4707 43772 9184 

PKS 0312-77 0.2230 0.1241 0.1999 37908 8651 

PKS 0405-12 0.5726 0.1241 0.5478 27208 7576 

HS0624-f6907 0.3700 0.1222 0.3464 61950 9184 

PG 0953+415 0.2390 0.1144 0.2163 24478 7747 

3C 249.1 0.3115 0.1222 0.2885 68776 9184 

PG1116-f215 0.1765 0.1144 0.1536 39836 8165/8097 

PG 1216+069 0.3313 0.1309 0.3078 69804 9184 

3C273 0.1580 0.1144 0.1417 18671 8017 

PG 1259+593 0.4778 0.1144 0.4533 95760 8695 

PKS 1302-102 0.2784 0.1183 0.2558 22119 8306 

PG 1444+407 0.2673 0.1222 0.2442 48624 9184 

3C 351.0 0.3719 0.1309 0.3483 73198 8015 

PHL1811 0.1917 0.1144 0.1690 33919 9418 

H 1821+643 0.2970 0.1144 0.2741 50932 8165 

Ton 28 0.3297 0.1231 0.3059 48401 9184 



^ The maximum redshift is defined for O VI absorbers at velocity separation 
> 5000 kms~^ from the background QSO, but the line search is conducted 
through the emission redshift of the QSO. 



orders were continuum normalized and co-added to form 
a single, stacked spectrum using our own software. The 
continuum was determined using a low-order polynomial 
fit to spectral regions that are free of strong absorption 
features. The combined echelle spectra of 16 QSOs have 
on average S/N > 7 per pixel over the majority of the 
spectral regions covered by the data. Together the total 
redshift survey path length covered by the echelle sample 
is A z 2.73 over z = [0.12, 0.5], comparable to the total 
re dshift path length at 2: < 0.15, Az — 2.2, carried out 
bv lDanforth fc Shulll (|2005h toward 31 sighthnes in FUSE 
spectra. Given the varying S/N versus wavelength along 
and between individual sightlines, we determine the A z 
versus the O VI 1031 absorption equivalent width thresh- 
old Wo in § 3.2. 

3. A STATISTICAL SAMPLE OF OVI ABSORPTION 
SYSTEMS 

3.1. The Search for VI Doublet Features 

To establish a statistically representative sample of 
O VI absorbers, we carry out a blind survey of O VI dou- 
blet features along each sightline that does not require 
known presence of hydrogen Lya or other ionic transi- 
tions. Before conducting the doublet search in each QSO 
spectrum, we first exclude spectral regions of known ab- 
sorption features from the Milky Way. Next, we search in 
the literature for known intervening absorbers, and fur- 
ther exclude spectral regions that contain these known 
features from the subsequent doublet search. For those 
sightlines not well studied, we performed this search our- 
selves. For the doublet search, we first smooth each 
stacked echelle spectrum using the Hanning function. 
Next, we step through the smoothed spectrum and search 
for deviant pixels that show negative fluctuation from 
unity by > 1.5 CTc, where CTc is the l-a error per smoothed 
pixel. We then determine whether a deviant pixel is a 
candidate absorption line based on a single gaussian- 
profile fit to the negative fluctuation. We restrict the 
width of the model profile to be no less than the spectral 
resolution clement and no larger than 300 kms^^. The 
last criterion is justified based on the line width distri- 



bution o f known O VI absorbers (|Heckman et al.l 120021 : 
iDanforth fc ShuU 2005). The deviant pixel is discarded 
from further consideration if the absorption equivalent 
width W based on the best-fit gaussian profile has < 2 a 
significance. 

Next, we treat each candidate absorption feature as 
the OVI 1031 line and determine the corresponding red- 
shift. Then we examine whether the spectrum is con- 
sistent with the presence of both members of the O VI 
doublet based on a simultaneous doublet profile fit at the 
expected locations of the redshifted OVI 1031 and OVI 
1037 lines. At this stage, we restrict both members to 
share the same line width and have a fiducial line ra- 
tio of 2 : 1, given the ratio of their respective oscillator 
strengths. The product of the doublet fitting procedure 
is a list of candidate O VI absorption doublets and their 
rest-frame absorption equivalent widths, 11/^^(1031) and 
VFr(1037). We visually examine every candidate and ac- 
cept doublets that satisfy both of the following two cri- 
teria. First, the OVI 1031 line has > 3cr significance. 
Second, the line ratio i?ovi = Wr(1031)/W^(1037) falls 
between 1 — au^ and 2 + 2 , where ct/Jq is the 

associated measurement uncertainty in i?ovi- Because 
the O VI 1031 transition is required to have at least 3-cr 
significance, the presence of both members together war- 
rants a higher significance for the detection of each O VI 
doublet. 

The entire sample of uncovered O VI absorbers include 
both foreground absorbers and absorbers that are as- 
sociated with the host of the background QSOs. In 
the following discussion, we define a "random sample" 
that consists of absorbers with |Ai;| > 5000 kms""'^ 
from the background QSOs and an "associated sam- 
ple" that are within |A-d| < 5000 kms"^ of the QSOs, 
where the ionization background is expected to be sig- 
nific antly different from the meta-galactic radiation field 
fe.g. iWevmann et al.|[T98lf ). The results of our search 
are summarized in Table 2. For each QSO line of sight, 
we report in columns (2) through (4) the total number 
of random absorbers, the minimum absorption equiv- 
alent width of the uncovered absorbers on that sight- 
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line, W^™, and the presence/ absence of absorbers in the 
vicinity of the background QSO. 

Figure 1 presents three examples of the O VI absorp- 
tion doublet uncovered in our blind search (panels a- 
c), together with an example of a previously reported 
absorber that is not recovered in our search (panel d). 
While our search routine has uncovered 27 absorbers to- 
ward 16 random lines of sight, four of 16 systems pre- 
viously reported toward seven lines of sight by other 
authors are not confirmed in our survey. We discuss 
these systems in more detail in § 5 and paper II. As 
a final step, we identify additional transitions associated 
with each fiducial O VI absorber. We perform a Voigt 
profile analysis to the O VI doublets, as well as Lya,/? 
(where available) using the VPFIT^ software package. In 
cases where absorbers have a complex structure, we at- 
tempt a multiple-component model profile. Note that we 
consider absorption components with velocity separation 
Av < 200 kms~^ as a single absorber in the statistical 
sample. 

3.2. Completeness of the Search 

To characterize the statistical properties of the OVI 
absorbers uncovered in our blind doublet search, it is nec- 
essary to first estimate the completeness of the survey. 
This is particularly important for the STIS data because 
(1) the S/N varies substantially across the entire wave- 
length range in each echcUe spectrum and (2) individual 
QSO spectra were obtained from different programs and 
span a wide range in their mean noise levels. We per- 
form a suite of Monte-Carlo simulations to quantify the 
completeness of our blind doublet search. A specific goal 
of the Monte-Carlo simulations is to estimate the total 
redshift path length A z covered by the survey for a given 
absorption equivalent width threshold in the OVI 1031 
transition, H^, , over which absorbers of Wr > can be 
recovered. Exploring the completeness versus allows 
us to maximize the siirvey efficiency, given the already 
small spectroscopic quasar sample. 

To carry out the simulations, we first mask re- 
gions along a single sightline k that are blocked by 
strong Galactic or extra-galactic absorption (e.g. Galac- 
tic 01 1302, or Galactic and extra-galactic Lya). Note 
that a strong absorption line at Xx will render us in- 
sensitive to doublet features in two regions in redshift 
space (Ax/1031.927 - 1.0 and Ax/1037.616 - 1.0), since 
it may mask either member of the O VI doublet. Next, 
we step through the redshift space defined by Zmin and 
zqso of the sightline defined in Table 1. At each redshift 

Zi, we adopt an intrinsic line width of 6 = 10 kms~ 
(corresponding to a gas temperature of T ~ 10^ K), gen- 
erate a model absorption line of varying strength Win, 
redshift the model to z^, and perturb the smooth model 
within the random noise characterized by the l-cr error 
array of the original echelle spectrum. Next, we apply 
the line search program and evaluate whether the syn- 
thetic absorption feature can be recovered. We repeat 
the procedure 1,000 times for each input line strength 
and determine the 3-cr detection threshold at Wq (zj), 
where the recovery rate of the synthetic absorption fea- 
tures is 99.7%. Finally, we repeat the procedure for all 
16 QSO spectra in the echelle sample. 

^ http://www.ast.cam.ac.uk/ rfc/ vpfit.html 



The result of the Monte-Carlo simulations is a sensitiv- 
ity curve (z) for every QSO sightline k that represents 
the 99.7% completeness of our blind OVI survey. These 
curves together allow us to calculate the total survey path 
length versus Wq for the entire sample of O VI absorbers 
from our search. For absorbers of Wr > Wq, they are 
expected to be detected at greater than 3-(T confidence 
level. We present in Figure 2 the collective Wo{z) curve 
from all 16 QSO sightlines together. At each redshift, the 
spread in Wq represents varying spectral quality in the 
echelle spectra of different QSOs. The plot shows that 
our blind search is sensitive to Wr ^ 30 mA absorbers 
at z « 0.2 — 0.5, while at z < 0.2 the majority of the 
spectra are sensitive to absorbers of Wr > 50 mA. 

To calculate the total survey path for out blind search, 
we first define a redshift path density g{Wr, z) at 

9{Wr, Zi) = J2 H{Z, - ^min) iJ(^max) H{Wr " 
fe 

(1) 

where H is the Heaviside step function, [z^j^, z^^^] are 
the minimum and maximum redshifts observed for the 
fcth QSO in Table 1, and the sum extends over all QSO 
sightlines in the echelle sample. The total redshift path 
length covered by the survey for a given equivalent width 
threshold Wr is therefore defined by 

A Z {Wr) = g{Wr) = j g{Wr, z) dz. (2) 

We evaluate equations (1) and (2) based on a red- 
shift interval of 5 z = 0.0001. The results are shown 
in panel (a) of Figure 3 for three different equivalent 
width thresholds, Wr = 30,50,80 mA. The plot illus- 
trates that, as stronger absorber are considered, the 
total survey path length increases; and that weaker 
absorbers require higher signal-to-noise data. Several 
strong dips are immediately apparent at redshifts z « 
0.18,0.21,0.22,0.26,0.29, which correspond to Galactic 
absorption from Ly a 1215, S II 1250, S II 1259, O 1 1302, 
C II 1334 respectively. 

In panel (b) of Figure 3, we show the redshift path 
along individual sightlines covered by our survey (solid 
black lines). Masked out regions due to Galactic ab- 
sorption features (e.g. the Lya 1215 transition) or other 
ionic transitions of known intervening absorbers appear 
as gaps in the path. Blue segments indicate the spectral 
regions that are sensitive to Wr = 30 mA. The tick- 
marks in panel (b) indicate the location of random O VI 
absorbers identified in our survey. We summarize the 
completeness of our O VI search in columns (5) and (6) 
of Table 2, where we list the total redshift path sensitive 
to Wr = 30 and 50 mA for each line of sight. 

4. STATISTICAL PROPERTIES OF THE O VI ABSORBERS 

We detect a total of 27 random O VI doublets along 16 
QSO lines of sight over a redshift range from zqvi = 0.12 
to zovi = 0.495 with a mean redshift of (z) = 0.25. 
Furthermore, we uncover associated absorbers toward 
nine QSOs at |Au| < 5000 kms^^ from the backgroimd 
QSOs. The random OVI absorbers form a statistically 
representative sample for constraining their global mean 
properties. The associated OVI absorbers form a uni- 
form sample for studying the gaseous halos of distant 
QSO hosts. Here we consider the statistical properties 
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Fig. 1. — Examples of O VI doublets. For each doublet, the A1031 member is plotted in the top panel and the associated A1037 line at the 
bottom. The corresponding 1-cr error array is shown in orange. We also present a best-fit Voigt profile of each absorber in blue. The spectra 
have been binned by a factor of 2 for display only. In panels (a) and (b), we present new OVI absorbers uncovered at z = 0.2028 and 
z = 0.3166 in our search routine, along the sightlines toward PKS0312— 77 and 3C351, respectively. In panel (c), we present a previously 
known absorber at z = 0.24532 toward H1821-f643 (Tripp ct al. 2000) that has also been recovered in our blind search. Finally, we show 
in panel (d) an absorber previously reported by these same authors at z = 0.1214 toward H 1821-1-643 based on additional FUSE spectra 
of the QSO. It is clear that the quality of available STIS spectra does not allow us to recover this absorber. Our random sample therefore 
does not include this system. 



B 




Fig. 2. — The 99.7% completeness limit in the absorption equivalent width, Wo, of the OVI 1031 transition versus redshift, for all 16 
sightlines in the STIS echelle sample. At each redshift, the spread in Wo represents varying spectral quality across different lines of sight. 
For absorbers at Zi with Wr > Woi^i), they are expected to be detected at greater than 3-cr confidence level. The plot shows that the 
STIS echelle sample as a whole is most sensitive to Wr ~ 30 mA absorbers at z fa 0.2 — 0.5. while at z < 0.2 the majority of the spectra 
are sensitive to absorbers of Wr > 50 mA. 
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Fig. 3. — Summary of the results from our blind search of OVI absorbers. In Panel (a), we present the total redshift path covered by 
the entire sample of 16 sighthnes for different absorption equivalent width thresholds Wt . Based on a series of Monte-Carlo simulations, 
we derive a total redshift path of A z = 1.66, 2.18, and 2.42 for Wr = 30, 50, 80 mA, respectively. In panel (b), we show for each individual 
QSO sightline the total redshift path allowed by the STIS echelle spectrum (black solid lines) and the redshift path that is sufficiently 
sensitive for detecting a Wr = 30 mA absorber at > 99.7% confidence level (blue lines). In panel (c), we present the number density of 
Wr > 30 mA OVI absorbers per sightline evaluated based on individual lines of sight, together with their associated errors. The errors are 
estimated according to Equation (4). We see a significant variation across the 16 sightlines that cannot be easily accounted for by random 
errors. 



TABLE 2 

Summary of the Blind O VI Survey and the Survey Limits 



QSO 

(1) 


jYrandom 

abs, 
(2) 


^min 

(mA) 
(3) 


Associated 
Absorbers 
(4) 


Az 

(W, = 50 m 
(5) 


{W > W*) 
A) (W. = 30 mA) 
(6) 


<W/dz 
{W* = 30 mA) 
(7) 


HE 0226-41101-2.3 


4 


43 


Y 


0.3055 


0.2593 


13 ±7 


PKS 0312-77 


1 


655 


N 


0.0375 


0.0014 


18 ± 18 


PKS 0405-124.5 


4 


30 


N 


0.3542 


0.2780 


12 ±6 


HS 0624+6907'' 


2 


27 


Y 


0.1681 


0.1458 


6±6 


PG 0953+415'^'* 


1 


121 


Y 


0.0698 


0.0434 


11 ± 11 


3C 249.1 


1 


72 


Y 


0.1206 


0.0826 


8±8 


PG 1116+215^.1° 


1 


76 


Y 


0.0307 


0.0053 


28 ± 28 


PG 1216+069 


1 


26 


N 


0.1339 


0.1047 


0± 10 


3C273" 


1 


22 


N 


0.0271 


0.0195 


0±51 


PGI259+59312 


2 


77 


N 


0.2851 


0.2680 


7±5 


PKS 1302-102 


2 


30 


N 


0.0992 


0.0596 


26 ± 19 


PG 1444+407 







Y 


0.0927 


0.0478 


0±21 


3C 351.0" 


1 


232 


Y 


0.1575 


0.1313 


6±6 


PHLI8III'' 


1 


64 


Y 


0.0388 


0.0133 


23 + 23 


H 1821+64315 


4 


30 


Y 


0.1226 


0.0999 


34 ± 17 


Ton 28 


1 


26 


Y 


0.1396 


0.0991 


+ 10 


total (16 sightlines) 


27 




10 


2.18 


1.66 


10 + 2 



Not e. — Sec also: fPILehner et al.l 1I2OO6'): (•2')ISavage et al.l ||2005|): (3') IGanguly et al.l (120061): (4) IChen &; P rochaskal 
I2OOOD: (5) jP rochaska ct al. (2004); (6) Aracil ot al. (2006); (7) Tripp & Savage (2000); (8) Sa vage et al.l ( 2002); (9) 
ISembach etal., ( , 2004); (1 0) Richter et al. ( 2 q06)_^_f 11) .Heap et al.i (i2002i) ; (12) iRichter et al., (,2004,) ; (Isi lYuan et al., (,200a) ; 
(14) IJenkins et al.l (|2003|) ; (15) iTripp'e't al.l | |2000() . 
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of the two samples separately. Studies of the physical 
nature of individual absorbers will be presented in a sep- 
arate paper (Thorn & Chen 2008; Paper II). 

4.1. The Incidence of Random O VI Absorbers 

The Monte-Carlo simulations presented in § 3.2 allow 
us to estimate the total redshift path length as a function 
of absorption equivalent width, from which we can then 
derive an accurate estimate of the number density of O VI 
absorbers per line of sight, using individual sightlines 
separately or the entire QSO echelle sample collectively. 
For each absorber of rest-frame OVI 1031 absorption 
equivalent width > W* identified in our search, the 
contribution to the total number density per line of sight 
dM/dz is l/g{Wl). Over a statistical sample, we can 
therefore derive for a given equivalent width limit VF,, 



dM{W > W^)/dz = ^ 



1 



X H{Wl.~W,) (3) 



where the sum extends over all systems found in a ran- 
dom sample and H is the Heaviside step function. Note 
that Az (W^) varies with W^. A stronger absorber may 
be detected in lower S/N spectra, allowing a larger sur- 
vey path length for a given QSO echelle sample. Equa- 
tion (3) is therefore not the total number of absorbers, 
^abs*^"™' divided by a fixed path length, (A z). The vari- 
ance of AM /dz is evaluated according to 



'dj^ /dz 



E 



[Az(W,0]2 



X H{Wi-W^). (4) 



We first evaluate AM / dz separately for individual lines 
of sight, adopting a common threshold for all QSO spec- 
tra. The results for W-j, — 30 mA are summarized in 
column (7) of Table 2, and presented in panel (c) of 
Fig. [3l Despite the relatively short path length per line 
of sight, it is immediately apparent from Fig. [3] that 
there is a substantial variation in AM jAz from one line- 
of-sight to another. In contrast to the line-of-sight to- 
ward H 1821-1-643, alo ng which four OV I absorbers are 
found over Az < 0.1 (jTripp et al.ll2000[ ). we note that 
the sightlines toward HS 0624-^6907, PG1259-h593, and 
3C 351.0 allow longer survey paths but exhibit fewer ab- 
sorbers (between 1 and 2). The large scatter in AM /Az 
between random sightlines underscores the necessity of a 
statistically significant sample of O VI absorbers for an 
unbiased estimate of their incidence. 

Combining all 16 QSO lines of sight together, we derive 
dM{W > 50 mA)/dz = 6 ± 3 over A z (M^* = 50 mA) = 
2.18, and dM{W > 30 mA)/dz = 9 ± 4 over A z {W^ ^ 
30 mA) = 1.66. In § 5 and panel (b) of Figure 5, we show 
that our number density measurements can be explained 
well by recent numerical simulations that do not require 
super starburst outflo w, but account for non-equilibrium 
ionization conditions (|Cen fc Fandl2006l ). 

4.2. The Cosmological Mass Density of Ionized Gas 

A particularly interesting quantity to determine, based 
on a statistically representative OVI absorber sample, 
is the cosmological mean mass density contained in the 
O^"*" ions, r2o^+- principle, we can then derive the 
cosmological mean mass density of warm baryons, fib, 
probed by the O VI transition according to 



toh 



f Zx (0/H)o mo 



(5) 



where / = iV(0'^+)/A^(0) is the fiducial ionization frac- 
tion of the OVI absorbing gas, Z = (O/H)/(O/H)0 rep- 
resents its mean metallicity, (0/H)q is the solar oxygen 
to hydrogen ratio by number, ^ is the mean atomic mass 
of the gas (/i — 1.3 for gas of solar composition), and 
mn and mo are the atomic weight of hydrogen and oxy- 
gen, respectively. In the following calculation, we adopt 
(0/H)o = 4.57 X 10"^ ( Asphmd et al. 2004). 

FoUowing lLanzetta et al.. (,1991.) , we calculate flo^+ ac- 
cording to 



CPc 
CPc 



E 



7V(0VI)/(iV(0VI))(iiV(0VI) (6) 
iV,(OVI) 



AX,, 



(7) 



where pc is the present-day critical density, f{N{0 VI)) is 
the frequency distribution function defined as the num- 
ber of O VI absorbers per column density interval per 
unit absorption path-length, Nt is the corresponding 
O VI absorption column density for absorber i with rest- 
frame absorption equivalent width and AXi is the 
co-moving absorption path length covered by the entire 
echelle sample at the absorption strength W^. The co- 
moving absorption path length A Xi for absorber i is 
related to the total redshift path length A z evaluated 
from Equation (2) at according to 

dX _ (l + z)^ 

~dJ ~ [Qm (1 + 2)3 - {Qm + - 1)(1 + z)2 + r!A]i/2 • 
The corresponding 1-a error is evaluated according to 



o5+ 



cpc 



E 



{Ax,y 



1/2 



(9) 



For each uncovered absorber of W^, we perform a Voigt 
profile fit to determine the corresponding Ni. 

We present the frequency distribution function deter- 
mined based on the random O VI absorber sample in 
panel (a) of Figure 5. The bin size is selected to contain 
roughly the same number of absorbers per bin. Error 
bars represent 1-a Poisson fluctuations. For the sam- 
ple of 23 random OVI absorbers with Wr > 30mA, we 
therefore derive ri(0^+) h = (1.7±0.3) x 10^"^. Assuming 
a Z = 1/10 Zq mean metallicity and / — 0.2 ionization 
fraction for the O VI gas, we further derive fib/i ~ 0.0014. 
Bearing in mind the large uncertainties in the metallicity 
and ionization fraction of the O^"*" bearing gas, our de- 
rived mean mass density of ionized gas is less than 10% 
of the total baryon mass density (e.g. Buries et al. 2001). 

4.3. Associated O VI Absorbers 

We have excluded spectral regions within |At;| — 
0—5000 kms^^ from the background QSOs for establish- 
ing the statistical sample of random foreground absorbers 
in the previous discussion, due to potential contamina- 
tion features from the QSO host environment. However, 
the blind OVI search has been carried out all the way 
through the emission redshifts of these QSOs. The OVI 
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Fig. 4. — Absorption profiles of associated OVI absorbers identified in our STIS echelle sample. 



absorbers found in the proximity of QSOs form a sample 
of associated absorbers that offers important constraints 
on the halo gas in QSO hosts (Chelouche et al. 2007, in 
prep). 

Ten of the 16 QSOs in our STIS echeUe sample exhibit 
strong OVI absorbers |Ai;| < 5000 kms~^ blue-shifted 
from their emission redshifts (indicated in Table 2). We 
present the absorption profiles of these OVI absorbers 
in Figure 4. The rest-frame absorption equivalent width 
of the 1031 member ranges from Wr = 60 niA (the ab- 
sorber at At; = -480 kms"^ from Ton 28) to = 3 A 
(the multi-component absorber at A u « —1600 kms^^ 
from 3C 351.0). The multi-component absorber found in 
the vicinity of 3C 351.0 has complex velocity profiles, ex- 
tending from Av ~ kms~^ to A?; « —2200 kms~^. 
It exhibits evidence of partial covering of the 0^+ ions 
and is understood to ori ginate in the o utflow gas from 
the background QSO ( Yuan et al.ll2002D . 

Excluding the complex sightline toward 3C351.0, the 
total redshift path length defined within \Av\ = — 5000 
kms^^ blue- ward of the emission redshifts of the remain- 
ing 15 QSOs in our echelle sample is Az = 0.33. At the 
equivalent width limit of our search W^, — 50 mA, we 
found nine O VI absorbers. We measure dJ\f/dz — 27±9 
for OVI absorbers of Wr(1031) > 50 mA in the vicinity 
of background QSOs. This is significantly higher than 
the incidence of random absorbers in the foreground. 



5. COMPARISONS WITH PREVIOUS SURVEYS 



Our search in the STIS echelle QSO sample has uncov- 
ered a total of 27 random O VI absorbers of Wr > 25 mA 
along 16 lines of sight, from which we derive the incidence 
of random OVI absorbers, dAf{W > 30 mA)/dz — 
10.4 ± 2.2 at (z) = 0.25. This is roughly half of what 
is mentioned in Tripp et al. (2006), who reported a to- 
tal of 44 OVI absorbers along 16 lines of sight and 
dJ\f{W > 30 mA) /dz = 23 ± 4 at z < 0.5. An in- 
dependent survey o f O VI ab s orber s in known Lya ab- 
sorbers by Danfort h fc Shulj ()2005f ) yields an absorber 
frequency of dJ\fiW > 30 mA) /dz = 17 ± 3 or dJ\f{W > 
50 mA)/dz = 9 ± 2 at z < 0.15. We present these differ- 
ent measurements in Figure 5, along with recent numeri- 
cal simulation results that include galactic wind feedback 
and non-equilibrium ionization conditions (|Cen fc Fand 
l2006f ). It is interesting to note that models with starburst 
feedback are not required by our measurements. 

The r esults of our blind searc h are consistent with 
those of 'Danforth & Jh^ (|2005[ ) from a FUSE survey 
at lower redshift that had been carried out in known 
Lya absorbers. The agreement is qualitatively consis- 
tent with the small fraction of OVI absorbers that do 
not have an associated Lya tra nsition. The discr epancy 
between our results and those of lTripp et al.l ()2006f ) based 
on a similar set of QSO spectral sample is most puzzling. 
In the following, we discuss possible reasons for this dis- 
crepancy. 

It is possible that the discrepancy is due to different 
QSO spectra covered by different surveys. Our random 
sample has been established based exclusively on avail- 
able STIS echelle data. The sightline along H 1821+643 
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has the highest incidence of O VI absorbers of aU 16 sight- 
Hnes in our echelle sample (See Table 2 and the discus- 
sion in § 4.1). In ad dition to th e four OVI absorbers 
originally reported in iTripp et al.l (l2000f) and confirmed 
in our blind search. IXripp et al.l ( 2001f ireported an addi- 
tional OVI absorber at z = 0.1212 in a targeted search 
in higher quality FUSE spectra. At z = 0.1212, the OVI 
A 1031 line is redshifted to Aobs ~ 1157 A, where the sen- 
sitivity of STIS cchcUe declines steeply. The lack of sen- 
sitivity at the blue end of the echelle sample is accounted 
for by the completeness study of our random sample dis- 
cussed in § 3.2. This z = 0.1212 absorber found in a 
targeted search after including additional FUSE spectra 
does not satisfy the criterion for establishing a statisti- 
cally unbiased sample, and has therefore been excluded 
in our random sample. 

In addition, there are OVI absorbers reported 
in the literature but not confirmed by our search. 
Along the same redshift path-length covered by our 
survey, 16 foreground O VI systems have been re- 
ported in the literature by other authors toward seven 
QSOs (HE 0226-4110, PKS 0405-123, PG 0953+415, 
PG 1116-^215, PG1259-I-593, H 1821-^643, and 3C273; 
See Table 2 for corresponding references), four are not 
confirmed in our survey. An example is shown in panel 
(d) of Figure 1, where we show the spectral region of an 
O VI dou blet along the s i ghtliii e toward PG 1259-1-593 re- 
ported bv lRichter et al.l i|2004D . The putative OVI 1031 
member is blended with the Si III 1206 transition of a 
Lya absorber aX z — 0.0461 (the feature at f = —30 
kms~^) and the GUI transition of a Lya absorbers at 
z = 0.2924 (the feature a.t v ^ +100 kms"^). Our 
search routine did not accept this absorber due to dis- 
crepant line strengths. Our random sample does not 
inclu de this possib l e abs orber. Along the same sight- 
line, iRichter et al.l (|2004[ ) also reported the presence of 
an OVI absorber at z = 0.3198, which is not confirmed 
by our routine due to discrepant absorption profiles and 
an inconsistent line ratio. 

Along th e sightline toward HE 0226-4110, 

iLehner et all (|2006l ) reported five OVI absorbers 
with Wr(1031) > 50 mA. We cannot confirm the pres- 
ence of their identification for a. z = 0.4266 absorber. In 
particular, the putative feature identified as the O VI 
1037 member is not present in our own reduction of 
the echelle spectra. Inspections of individual exposures 
showed that the broad feature are predominant in two 
of the 12 exposures, but the discrepancy may be due to 
different treatment in the hot pixels in echelle spectra 
(T. Tripp, private communication). Additional spectra 
obtained using the Gosmic Origin Spectrograph (GOS) 
or STIS will be valuable for confirming the presence of 
the doublet. Gonsequently, this absorber is not included 
in our random sample. 

6. DISCUSSION 

6.1. The Ionization State of the O VI Gas 

Various studies of individual O VI absorbers that con- 
sider both kinematic signatures of the absorption pro- 
files and abundance ratios between different ions have 
concluded that a large fraction of these absorbers can 
be explained by a simple photo-ionizat i on model (e.g. 
iTripp et al.ll2000t iProchaska et all 12004 ISembach etall 




Fig. 5. — (a) The frequency distribution function of OVI ab- 
sorbers from our random sample. The bin size is chosen for each 
bin to contain roughly the same number of absorbers. Error bars 
represents l-cr Poisson fluctuations, (b) Observed incidence of O VI 
absorbers from differ ent studies. Star points are from this work; 
open circles are from 'Danforth & Shull (20M) ; and open squares 
are from Savage et al. (2002); Prochaska et al. (2004); Tripp et ^ 
(120061 ): iCooksev et al.i (i2008i . Note that we have slightly offset 
the data points in the horizontal direction for clarification. The 
curves are recent WHIM predictions based on numerical simula- 
tions (Ceil & Fang 2006) that account for non-equilibrium ioniza- 
tion conditions. The top curve further includes Galactic wind feed- 
back. 



l2004HLehner et al.ll2006l ). These results are in contradic- 
tion to the expectation of a WHIM origin, under which 
the bulk of the 0^+ ions originate in the shocks created 
as the gas accretes onto the large filamentary structures. 
However, it is often difficult to rule out a coUisional ion- 
ization scenario for these absorbers. Additional uncer- 
tainties arise, if OVI absorbers arise in post-shock gas 
that is undergoing radiative cooling and a simple ion- 
ization equilibrium assumption does not apply (Gnat & 
Sternberg 2007). 

The comparison of the absorption column densities 
in Figure 6 between Lya and O VI transitions for the 
OVI selected sample (solid points) demonstrates that 
the A^(H I) associated with these OVI absorbers spans 
nearly five orders of magnitude. In addition, stronger 
O VI transitions have on average stronger associated Lya 
transitions. The null hypothesis in which Lya and OVI 
column densities are randomly distributed among each 
other can be ruled out at > 97% confidence level based 
on a generalized Kendall test that considers the non- 
detection of Lya for the O VI absorber at z = 0.32639 to- 
ward HE 0226-4110. We note that the STIS echelle sam- 
ple is sensitive to O VI absorbers of Wr > 80 mA over 
the entire redshift path-length allowed by the data, and 
O VI absorbers of Wr > 50 mA over « 80% of the sight- 
lines. The echelle spectra therefore allow us to uncover 
OVI absorbers of log 7V(0 VI) > 13.5 over the major- 
ity of the sightlines. The correlation is therefore unlikely 
due to a selection bias. Including measurements from 
known O VI absorbers at low redshift (open squares and 
crosses), we find that the correlation remains at > 80% 
significance level. 

The broad range of A^(H I) spanned by these O VI ab- 
sorbers also makes it difficult to apply a simple coUi- 
sional ionization equilibrium model for explaining their 
origin. Furthermore, the correlation between N{OVl) 
and N{}i I) indicates that the OVI absorbers appears to 
track closely HI gas in ov erdense regions, while t he op- 
posite does not apply fe.g. lDanforth fc Shullll2005t ) . This 
trend can be understood under a simple photo-ionization 
scenario, and the slope in A^(0 VI) versus N(H I) implies 
a declining ionization parameter (defined as the number 
of ionizing photons per gas particle) with increasing gas 
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column density (see Figure 31 of Paper II). 

6.2. O VI Absorbers without Associated HI 

A unique feature of our bhnd search is that the find- 
ing of O VI doublet does not rely on known presence of 
Lya transitions. The sample therefore allows us to con- 
strain the fraction of OVl absorbers that do not show 
significant traces of HI. These absorbers may represent 
an extreme case for outfiow gas that has been ejected 
into underdense region s by strong galactic winds (e.g. 
iKawata fc Rauchll2007l ). 

In Figure 6, we find that 26 of the 27 O VI absorbers in 
the random sample show associated HI absorption fea- 
tures. The only one possible Lya-free system is found 
to have iV(0 VI) = 13.6 ± 0.1 at zqvi = 0.3264 toward 
HE 0226— 4110 (the arrow in Figure 6). The comparison 
indicates that < 5% of OVI absorbers originate in un- 
derdense regions of Ar(H I) < lO^'^ cni^^ that are heavily 
metal-enriched by starburst outflows 

6.3. Large-scale Variation in Absorber Number Density 

Cosmological simulations also show that OVI ab- 
sorbers originating in the WHIM reside in regions of 
over-density 6 = p/ (p) = 10 — 100 at z ^ (jPave et al.l 
[200llCen et"all[2?)0ll : iFang fc Brvanl[200l . This is sim- 
ilar to what is found in simulation studies of Lya forest 
absorbers, according to which by z = the majority 
of Lya absorbers are associated with 6 — 10 — 100. It 
is therefore interesting to examine the relation between 
Lya and O VI absorbers. 

Of the 16 QSO sightlines in our echelle sample, 
seven have been studied extensively f or the line-of- 
sight properties of Lya absorbers (ISembach et all 
[2003 iWilliger et aP [20061 : iLehner et al.l 120061) . We 
present in Figure 7 a comparison of the observed 
number densities between Lya absorbers and O VI 
absorbers measured along individual sightlines. These 
include H 1821-^643, HE 0226-4110, HS 0624-1-6907, 
PG 0953-^415, PG 1116-1-215, PG 1259-^593, and 
PKS 0405— 123. The statisti cs of Lya absorber s are 
taken from the compilation of 'Lehn er et al.l (|2006'). We 
have included the corrected Lya absorber catalog for 
the PKS 0405— 123 sightlinc. The error bars represent 
1-cr uncertainties. 

Considering all seven sightlines together, we find a 
mean number density of {dM/dz) = 66 for the Lya 
absorbers with an r.m.s. dispersion of 15. In contrast, 
we find {dJV/dz) = 12 for the OVI absorbers with an 
r.m.s. dispersion of 10. Examining the number density 
variation from one sightline to another, we find that the 
line of sight toward H 1821-1-643 with the highest inci- 
dence of O VI absorbers has the lowest incidence of Lya 
absorbers. Figure 7 shows a possible inverse correlation 
between the number density of H I absorbers and the 
number density of OVI absorbers. While Lya absorbers 
appear to be relatively more abundant than O VI along 
each sightline, the number density of O VI absorbers do 
not scale proportionally with the number density of Lya 
absorbers along individual lines of sight. 

This anti-correlation between [dJ\f/dz]ovi and 
[dj\f / dz]-^y.^ may be understood, if those sightlines with 

a higher incidence of O VI absorbers probe a larger 
fraction of more highly ionized regions. However, we 



note the result from § 6.1, where we find nearly every 
one of the O VI absorbers in our statistical sample shows 
an associated Lya transition. On the other hand, if 
O VI absorbers originate preferentially in underdense 
regions, then a higher incidence of O VI implies a 
substantial fraction of the sightline probe underdense 
regions and hence a lower incidence of Lya absorbers. 
This is expected if heavy elements from starburst 
outfl ows are ejected primari ly into underdense regions 
(e.g. IKawata fc RauchI r2007f ). Further insights can be 
obtained from a detailed examination of the large-scale 
galaxy environment of the OVI absorbers. A larger 
sample of QSO echelle spectra will also be valuable for 
a more quantitative study to confirm this result. 

7. SUMMARY 

We have conducted a systematic survey of intervening 
O VI absorbers in available echelle spectra obtained using 
HST/STIS with the E140M grating. We have collected 
echelle spectra of 16 distant QSOs from the HST data 
archive. These UV bright quasars at xqso = 0.17 — 0.57 
have been selected for HST/STIS observations by dif- 
ferent authors for different studies, and represent a ran- 
dom, flux-limited sample of QSO echelle spectra for an 
unbiased search of intervening absorbers. The primary 
objectives of our blind OVI survey are: (1) to obtain 
a statistically representative sample of O VI absorbers 
along random lines of sight for an accurate measure- 
ment of the cosmological mass density of warm-hot gas 
at low redshift; and (2) to study the incidence of OVI 
absorbers along individual sightlines for quantifying the 
spatial variation of OVI absorbing gas. 

Our search in the STIS echelle QSO sample has un- 
covered a total of 27 random {\Av\ > 5000 kms~^ 
from the background QSO) O VI absorbers with rest- 
frame absorption equivalent width of the 1031 mem- 
ber Wr ^ 25 mA along 16 lines of sight. Ten of these 
QSOs also exhibit strong OVI absorbers at \Av\ < 5000 
kms^ blue-shifted from the QSO emission redshifts. 
Our O VI survey differs from previous studies in two fun- 
damental ways. First, we do not require prior knowledge 
of the presence of other transitions such as Lya. The 
echelle resolution allows us to identify the O VI absorp- 
tion doublet based on their common line centroid and 
known flux ratio. Second, the survey completeness is 
well understood and characterized based on a series of 
Monte-Carlo simulations. We obtain an equivalent width 
threshold curve for every QSO sightline. Together, these 
curves allow us to estimate the total redshift survey path 
as a function of equivalent width threshold, thereby max- 
imizing the survey efficiency at the limit allowed by the 
echelle spectra. 

Based on the results of the Monte-Carlo simulations, 
we derive a total redshift path for random absorbers of 
Az = 1.66,2.18, and 2.42 for Wr ^ 30,50 and 80 mA, 
respectively. The total path-length surveyed in our sam- 
ple is comparable to the FUSE survey carried out at 
z < 0.15 by Danforth fc ShuU (2005), and is the largest 
at z = 0.12 — 0.5. The estimated total redshift path leads 
to a number density of dJ\f{W > 50 mA) /dz = 6.7 ±1.7 
and dAf(W > 30 mk)/dz = 10.4±2.2. Recent numerical 
simulations that allow non-equilibrium ionization condi- 
tions explain the observations very well. In contrast, 
we also measure dN /dz = 27 ± 9 for OVI absorbers of 
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Fig. 6. — Comparison between A'^(0 VI) and A'^(H I) for the 27 O VI absorbers identified in our bUnd search (solid points. All absorbers but 
one have an associated Lyg abso rption line. We have also inclu ded published absorbers at low reds hifts. Open squares a re z < 0.12 absorbers 
presented in 'Savage et al.' (200^) for PG 0953+415, in ScmbadTct al. (2004) for PG 1116+215, in lRichter eTatl l l200i) fo PG 1259+593, in 
[Prochaska ct al. (2004) for PKS 0405-123, and in .Cookscv ct al. (2008.) for PKS 1302-102. Crosses indicate measurements from a FUSE 
survey bv lDaniforth Shulll I I2005I) . 
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Fig. 7. — Comparison of the line density variation from one line of sight to another for hydrogen Lya absorbers and O VI absorbers. The 
number density of Lyo absorbers measured for individual sightlines are taken from Lehner et al. (2007) for Lya absorbers of rest-frame 
equivalent width Wr'(Lya) > 90 mA and the number density of OVI absorbers measured for individual sightlines are from the blind search 
for O VI absorbers of rest-frame equivalent width Wr{0 VI) > 30 mA presented in this paper. We find a possible inverse correlation between 
the number density of H I absorbers and the number density of O VI absorbers. 
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VFr(1031) > 50 mA in the vicinity of background QSOs. 
For the sample of 27 random O VI absorbers, we further 
derive a mean cosmological mass density of the O VI ions 
n(0^+) h = (1.7 ± 0.3) X 10"'^ that leads to a mean cos- 
mological mass density of warm baryons f2b^ ~ 0.0014. 

Using the statistical sample of O VI absorbers, we show 
that < 5% of O VI absorbers originate in underdense re- 
gions that do not show significant trace of HI. In ad- 
dition, a comparison of the absorption column densities 
between Lya and OVI transitions for the OVI selected 
sample shows that the A'^(H I) associated with these O VI 
absorbers spans nearly five orders of magnitude, making 
it difficult to apply a simple coUisional ionization equi- 
librium model for explaining their origin. Furthermore, 
stronger OVI transitions have on average stronger as- 
sociated Lya transitions, indicating that the OVI ab- 
sorbers appear to track closely H I gas in ovcrdcnsc re- 
gions, while the opposite does not apply. At the same 
time, comparisons of individual line-of-sight properties 
show a moderate inverse correlation between the num- 
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